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The importance of Hund's rule coupling for the stabiliza- 
tion of itinerant ferromagnetism is investigated within a two- 
band Hubbard model. The magnetic phase diagram is calcu- 
lated by finite-temperature quantum Monte Carlo simulations 
within the dynamical mean-field theory. Ferromagnetism is 
found in a broad range of electron fillings whereas antifer- 
romagnetism exists only near half filling. The possibility of 
orbital ordering at quarter filling is also analyzed. 



I. INTRODUCTION 

Investigations of the microscopic mechanisms responsible 
for the stability of metallic ferromagnetism have recently re- 
ceived renewed attention. In particular, it was finally estab- 
lished that itinerant ferromagnetism is indeed stable in the 
one-band Hubbard model at intermediate on-site interactions, 
an important condition being a properly tuned kinetic energy 
with a pronounced peak in the density of states (DOS) near 
the band edge [Q-^. For real ferromagnets, e.g., the transition 
metals Fe, Co, and Ni, quite a different microscopic mecha- 
nism for ferromagnetism, based on Hund's rule coupling in the 
presence of orbital degeneracy, is also expected to be relevant. 
It was first suggested by Slater and then stressed by van 
Vleck [0 that the intra-atomic exchange leading to "Hund's 
rule atomic magnetism" might be responsible for bulk ferro- 
magnetism, i.e., the hopping of electrons or holes might lead 
to a bulk ordering of preformed atomic moments. 

Concrete calculations based on microscopic models began 
with Roth who considered a two-band extension of the 
Hubbard model 

+ ^ {Vo — 5^„,Fo)ni„uni^' ct' 



(1) 



where v denotes the two (or more) orbitals, a the spin, and i, j 
the lattice sites. In this model all interactions are purely local, 
i.e., occur only on a single site. Apart from a Hubbard inter- 
action U for electrons of opposite spin on the same orbital, 
there is also a density-density interaction Vq between electrons 
on different orbitals, as well as an intra-atomic exchange in- 
teraction _Fo which is separated into its density-density and 
spin-flip contribution (the spin-flip term is neglected in Sec. 
III). An "on-site pair hopping" -term of the same size is not 
considered since it requires an empty and a doubly-occupied 
orbital to take effect. Such conflgurations are strongly sup- 
pressed by the Hubbard interaction U and the kinetic energy. 



The kinetic energy describes the hopping of electrons of a 
given spin between identical orbitals on different sites. 

It is instructive to consider the two-site model with one 
electron per atom, i.e., at quarter filling. In the strong- 
coupling regime the ground state is a spin triplet and orbital 
singlet, i.e., on the two atoms different orbitals are occupied. 
Within second order perturbation theory the energy of this 
state is readily calculated as —At^/iVo — Fq). This connec- 
tion between staggered orbital ordering and ferromagnetism 
was investigated already in 1966 by Roth ||l^: applying the 
random phase (or Hartree-Fock) approximation to Km she ob- 
served that, at r = and quarter filling, a ferromagnetic 
state is unstable against an additional staggered orbital or- 
dering for Vo — (7/2 — Fo/2 > 0. For decreasing temperatures 
the Hartree-Fock approximation predicts first a phase transi- 
tion from paramagnetic to ferromagnetic ordering, and at a 
lower temperature a second transition to a phase with ferro- 
magnetic and orbital ordering. We note that the phenomenon 
of staggered orbital ordering may be a characteristic feature 
of any effective two-band model with orbital degeneracy (as 
realized, for example, in eg-bands). 

A major step towards the understanding of the physics of 
the two-band model was the derivation of an effective strong- 
coupling Hamiltonian at quarter filling by Kugel' and Khom- 
skii |l^, and by Cyrot and Lyon-Caen |l^ who included the 
effect of on-site pair hopping. This effective Hamiltonian has 
coupled spin and orbital (pseudo spin) degrees of freedom. 
Within a self-consistent-field approximation |l^, where or- 
bital and spin degrees of freedom are decoupled, and also 
within the molecular field theory |l4|, where correlations be- 
tween different lattice sites are neglected, the effective Hamil- 
tonian shows an insulating ferromagnetic ground state with 
staggered orbital ordering. Contrary to the weak-coupling 
Hartree-Fock approximation this strong-coupling approach 
predicts that orbital ordering occurs first when the temper- 
ature is decreased jl4|. The ferromagnetic ground state at 
quarter filling and sufficiently strong Coulomb interactions 
was confirmed by exact diagonalization studies of finite sys- 
tems in one dimension [|l5|-|l7|. 

Off quarter or half filling the high degeneracy in the atomic 
limit makes a perturbational analysis essentially impossible 
[Q. Nevertheless, for infinite Coulomb interaction ?7 = oo 
and one dimension, Miiller-Hartmann |^ proved that for 
Fo > the ground state of (|l|) is ferromagnetic for fillings 
< n < 2. Recently exact diagonalization studies in one 
dimension were performed by Hirsch The results de- 

pend sensitively on the boundary conditions and the number 
of lattice sites: below quarter filling and for six lattice sites 
ferromagnetism is found for the system with open boundary 
conditions, while it is absent when periodic boundary con- 
ditions are used. Insight was also obtained by several ap- 
proximative treatments, in particular the Hartree-Fock the- 
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ory [ [L4yi8| ] , a generalized Hartree-Fock approach [gO|, more 
complex variational wave functions j2^, the local approach 
[ p2| , and the Gutzwiller approximation [^. They all find 
ferromagnetism to be stabilized by Hund's rule coupling at 
intermediate to strong Coulomb interactions. Clearly, in this 
regime a proper treatment of correlation effects and the dy- 
namics of the quantum mechanical many-body problem is es- 
sential. In the last few years the dynamical mean-field theory 
(DMFT) has turned out to provide a reliable framework and 
powerful method for the investigation of such types of prob- 
lems. Hence we use it in this paper to study the stability of 
the ferromagnetic phase at and off quarter filling for inter- 
mediate values of the Coulomb interaction, and to determine 
transition temperatures. 

The paper is structured as follows. In Sec. ^the quantum 
Monte Carlo (QMC) algorithm to solve the DMFT equations 
is introduced. The magnetic phase diagram is presented in 



Sec. 



[II A 



in Sec. |IIIB 



and the possibility of orbital ordering is discussed 



II. QUANTUM MONTE CARLO SOLUTION OF 
THE DYNAMICAL MEAN-FIELD EQUATIONS 



The DMFT P4H30| approximates the lattice model by a 
single-site problem of electrons in an effective medium (mean- 
field) that may be described by a frequency dependent, i.e., 
dynamical, self-energy Syo-(w). The latter has to be deter- 
mined self-consistently via a k-integrated Dyson equation that 
reads 



Guan — I ds — 



(2) 



with Matsubara frequencies uj„, Green function G^o-n ~ 
GuaiiuJn), and the DOS of the non- interacting electrons 
N°{e). In the case of the multi-band Hubbard model (^) the 
single-site problem takes the form 

G..„ = -J J V[mr]i>..ur..ne^^^'^'-'^~"^ (3) 

where ip and ip* are Grassmann variables, = + E, 
and Al^pj^p* ,Q~^] denotes the single-site action including all 
local interactions 
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As for the one-band model ||^,||| Eqn. (|) is equivalent to 
a (degenerate) Anderson impurity model and can thus be 



treated by standard techniques. In the present paper the 
QMC algorithm of Hirsch and Fye jsi] will be used. In a first 
step each one of the four terms of Eqn. (^ are decomposed via 
the Trotter-Suzuki formula, and imaginary time is discretized 
(r = I At for I = 1 ... A and At — (3/A). In the second step 
the interaction terms are decoupled to obtain a quadratic ac- 
tion. For the density-density interactions this is achieved as 
usual by the Hubbard-Stratonovich transformation 



exp{ — (Vo - Fo){:4}l<rli>ual - ijt'altp^'al) } = (5) 



where cosh(A) = exp(Ar(Vb — Fo)/2). The same decoupling 
holds for the remaining Vo and the U terms. For A'^ or- 
bitals this yields altogether {2N — N)A auxiliary fields for 
the density-density interactions. 

By contrast, there does not exist a standard decoupling 
scheme for the last term of Eqn. (^ describing a spin-flip. Re- 
cently Motome and Imada j32j proposed a decoupling scheme 
for this term: it avoids the minus-sign problem in the sym- 
metric case at half filling, but leads to a phase problem off half 
filling due to the use of a complex auxiliary field. Instead, we 
tried the following Hubbard-Stratonovich decoupling: 



exp{ATFoil;t„iili„^„iil!*,_^iip^,^i) = 



(6) 



with /i — \/ At Fo ■ Unfortunately, this transformation was 
found to lead to a sign-problem, too, even within the DMFT 
where the sign-problem is absent in the single-band case. 
For this reason we neglect the spin-flip term in the follow- 
ing. Hund's rule coupling is thus restricted to the direction 
of the quantization axis, i.e., the 2-axis, implying the break- 
ing of the SU(2) spin rotation symmetry. The restriction to 
an Ising-type Hund's rule coupling is commonly used in the 
investigations of the two-band Hubbard model |l2[[l^ . |l9|j2c[| . 
This restriction has no effect on the critical temperatures in 
the limits of weak and strong coupling, i.e., in the Hartree- 
Fock approximation and the Weiss mean-field theory of the 
corresponding spin model, respectively. Therefore we expect 
that the spin-flip term has no strong effect on the critical tem- 
perature at intermediate coupling. Note, that while the influ- 
ence of the spin-fiip term on critical temperatures is probably 
small the excitation spectrum at T=0 depends sensitively on 
this term. It shows a spin gap or not, respectively. 

Based on the above Hubbard-Stratonovich decoupling the 
Monte Carlo method is employed to sample the auxiliary 
spin configurations and thus to calculate Green functions and 
susceptibilities. The multi-band algorithm is a generaliza- 
tion of the one-band algorithm. A similar generalization was 
employed by Rozenberg |Q in the investigation of metal- 
insulator transitions in the two-band Hubbard model without 
Hund's rule coupling. 

Phase boundaries are determined by a Curie- Weiss fit of 
the corresponding susceptibility in the homogeneous phase. 
Since the corresponding critical temperatures still depend on 
the unphysical time discretization parameter At a second or- 
der fit to At = was performed from at least six values of 
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At G [0.075,0.5]. Besides the statistical error of the QMC 
simulation (propagated via T- and Ar-fit), there exists an 
additional systematic error due to higher order contributions 
in the T and At fits. In particular the At dependence of the 



III 



Curie temperature was considerable for the data of Sec 
Therefore we estimated this systematic error by comparing 
the Curie temperature obtained from all At values to that 
calculated without the At = 0.5 value. The individual and 
mean difference between these two fits was within the sta- 
tistical error of the Monte Carlo data. This analysis implies 
that the systematic error is smaller than the statistical error 
shown in the figures. 



III. RESULTS AND DISCUSSION 
A. Ferromagnetism 

In the present paper we investigate the Hund's rule cou- 
pling in the presence of orbital degeneracy as a possible ori- 
gin of itinerant ferromagnetism. This microscopic mechanism 
should be distinguished from the one found to be important 
for the single-band Hubbard model, which is based on an 
asymmetric DOS [Q-^. Therefore we employ a symmetric 
semielliptic DOS N°{e) = ^(2t*)2 -e2/(27rt*2) (in the fol- 
lowing t* = 1 will set our energy scale). For this DOS and 
a symmetric Gaussian DOS no ferromagnetism was found in 
the single-band Hubbard model up to a Coulomb interaction 
of f/ = 30 |,|. 

In fact, even the two-band model without Hund's rule cou- 
pling does not indicate ferromagnetism at U — 9, Vo = 5, 
and n = 1.25 (see Fig. |l]). However, already a small Hund's 
rule coupling, Fp = 1.61 ± 0.15, is sufficient to stabilize ferro- 
magnetism |34|. This value is considerably smaller than that 
obtained by Kuei and Scalettar |l6j: in one dimension (six 
sites) at n = 1 with U = Vo + Fq no ferromagnetism was 
found below Fo ~ 8 (in units rescaled to obtain a bandwidth 
of four). This value of Fo is, however, comparable to the 
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FIG. 1. Curie temperature vs. Hund's rule coupling Fo for 
the two-band Hubbard model at 17 = 9, Vb = 5, and n = 1.25. 
Ferromagnetism is seen to be stabilized by Hund's rule cou- 
pling. 




FIG. 2. Magnetic T — n phase diagram for U = 9, 
Vo = 5, and Fo = 4 with paramagnetic (PM), ferromagnetic 
(FM), and antiferromagnetic phase (AFM). Ferromagnetism 
is found in a broad range of fillings. 




FIG. 3. Same as Fig. g but calculated within the 
Hartree-Fock approximation (solid line: Curie temperature, 
dashed line: Neel temperature, dotted line: first order transi- 
tion between FM and AFM). 

results of Hirsch [0: in one dimension at n = 1 and n = 
1.5 (six and four sites, respectively) with U = Vo -\- 2_Fo no 
ferromagnetism was found below Fo ~ 1.5 and 2, respectively. 
Note, that the relation 17 = Vb + Fo makes the Hamiltonian 
(0) form-invariant w.r.t. orbital rotations, while this holds 
for U = Vo + 2Fo if an additional pair hopping term is added. 
The Hund's rule coupling Fq must be smaller than the density- 
density interaction Vo since otherwise an unphysical attractive 
interaction between two electrons on the same site (in different 
orbitals with the same spin) exists. 

The magnetic T~n phase diagram is presented in Fig. ^for 
17 = 9, Vb = 5, and a relatively strong Hund's rule coupling 
Fo = 4. Ferromagnetism is found in a broad range of elec- 
tron densities n. Note, that the phase-diagram is symmetric 
around n = 2 due to particle-hole symmetry. The special case 
of quarter filling n — 1 does not mark a pronounced point in 
the phase diagram. In particular the maximal Curie temper- 
ature is found above quarter filling. Typical Curie tempera- 
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tures are about 0.1, which for a bandwidth of 4eV {t* = leV) 
corresponds to O.leV, i.e., about lOOOK. Near half filling the 
antiferromagnetic Heisenberg exchange suppresses the ferro- 
magnetic order and a narrow antiferromagnetic phase with 
the usual AB sublattice structure develops. Fig. ^ suggests 
that Curie and Neel temperature cross at a finite tempera- 
ture and not at T = 0. This is confirmed by the observation 
that for n = 1.85 the Curie temperature extrapolated from 
the paramagnetic phase lies only slightly below the Neel tem- 
perature. 

The same T — n phase diagram calculated within the 
Hartree-Fock approximation is shown in Fig. ^. The Hartree- 
Fock approximation fails to describe the suppression of the 
magnetic order at the crossover from ferro- to antiferromag- 
netism. Furthermore, the magnetic phases are overestimated, 
i.e., the critical temperatures are more than an order of mag- 
nitude too large and both magnetic phases continue to ex- 
tremely small values of n. Within the Hartree-Fock approx- 
imation a first order phase transition between antiferromag- 
netic and ferromagnetic phase occurs. It is not clear at present 
whether this is also true within the DMFT. 




FIG. 4. Critical temperature for the instability against a 
state with mixed ferromagnetic and orbital ordering (control 
parameter 5) as obtained from [x(7c,<5)]~^ = 0. For 17 = 9, 
Vb = 5, -Fo = 4, and n — 1 the highest temperature is found 
at 5 = 0, indicating a transition into a purely ferromagnetic 
state. 



B. Orbital ordering 

As was pointed out in the introduction, at quarter filling the 
Weiss mean-field theory for the effective Kugel' and Khomskii 
Hamiltonian [ p"3|Jl^ predicts an instability of the paramag- 
netic phase against staggered orbital ordering when the tem- 
perature is decreased, while a transition to pure ferromagnetic 
order is suggested by the Hartree-Fock approximation. The 
DMFT, containing both approximations as limits at strong 
and weak coupling, respectively, is well suited to clarify this 
contradiction. To investigate a possible transition from the 
phase without long-range order to a phase with mixed ferro- 
magnetic and staggered orbital ordering we introduce a pa- 
rameter 5 that allows one to investigate the instability w.r.t. a 
general mixed ordering. For every 6 the corresponding sus- 
ceptibility is defined through a field h, which modifies the 
grand potential Q, by the term 



1 d^Q, 
L'dh?' 



{(l-5)a + 5(-l)V} 



(7) 
(8) 



Here L denotes the number of lattice sites. For 5 = one 
obtains the (para-) magnetic susceptibility and for 5=1 the 
susceptibility for orbital ordering. Fig. ^shows that for U = 9, 
Vo = 5, and Fo = 4 the critical temperature is maximal for 
5 = 0, i.e., the paramagnet is unstable against pure ferro- 
magnetic order. Thus, the phase diagram Fig. ^ need not be 
modified by phase transitions from paramagnetic to orbital 
ordering. Whether the ferromagnetic phase becomes unsta- 
ble against orbital ordering at even lower temperatures can- 
not be answered by the method employed here, i.e., by the 
calculation of susceptibilities within the paramagnetic phase. 

The scenario of Fig. ^qualitatively agrees with the Hartree- 
Fock approximation, which even at T = does not predict 
orbital ordering since Vo — U/2 — Fo/2 < 0. The disagreement 
with Weiss mean-field theory can be explained by the fact 




0.4 0.6 0.8 

5 

FIG. 5. Same as in Fig. § but for [/ = 8, l^o = 6, i^o = 2, 
and n = 1. Here, the maximal temperature is found at 5 = 1, 
i.e., an instability against staggered orbital ordering occurs. 

that the second order perturbation theory leading to the 
effective strong-coupling Hamiltonian is controlled in t'^ /{Vo~ 
Fo) which is of 0(1) here. 

The considerations above suggest that orbital ordering may 
occur if the intra-atomic exchange Fq becomes smaller. In- 
deed, for U — 8, Vo = 6, and Fo — 2 a phase transition to pure 
orbital ordering (5=1) occurs (see Fig. Fig. ^ shows the 
T — n phase diagram for these parameters, including orbital 
ordering near quarter filling, antiferromagnetism near half fill- 
ing and ferromagnetism in between. We should mention, that 
although at n = 0.9 the inverse orbital ordering susceptibility 
decreases with decreasing temperatures in a Curie- Weiss like 
behavior, suggesting a transition at T « 0.04, it decreases 
again at lower temperatures. Hence orbital ordering does not 
take place at n = 0.9. 

Even at n = 1.2, i.e., the crossing point of the orbital order- 
ing and ferromagnetic phase boundary in Fig. no instability 



4 



against a mixed ferromagnetic and orbital ordering is found. 
Therefore we conclude that the two-band Hubbard model 
with Hund's rule coupling Fq shows an instability towards 
either pure ferromagnetic or pure orbital ordering. However, 
since these phases do not exclude each other a phase with 
mixed order may appear at even lower temperatures, as is 
predicted by Hartree-Fock and Weiss mean-field theory. 




FIG. 6. T—n phase diagram for U = 8, Vo = 6, and Fo — 2. 
In addition to the phases in Fig. ^ an orbital ordering (OO) 
phase is found near quarter filling. 



Conclusion 

We showed that even for a symmetric DOS the Hund's rule 
coupling provides an effective microscopic mechanism for the 
stabilization of ferromagnetism in a broad range of electron 
densities off half filling. This mechanism takes effect at in- 
termediate to strong values of the Coulomb interaction and 
therefore requires a proper treatment of the quantum mechan- 
ical correlations. It is different from the mechanism based on 
an asymmetric DOS, which leads to ferromagnetism even in 
the single-band Hubbard model. The question which one of 
these mechanisms is the main driving force for ferromagnetism 
in Fe, Co, and Ni remains open, if it can be explained by a 
single mechanism at all. 
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